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ABSTRACT 

' We present an analysis of the properties of the hot interstellar medium (ISM) in the merging 

pair of galaxies known as The Antennae (NGC 4038/39), performed using the deep, coadded 
^411 ks Chandra ACIS-S data set. These deep X-ray observations and Chandra's high angular 
resolution allow us to investigate the properties of the hot ISM with unprecedented spatial and 
spectral resolution. Through a spatially resolved spectral analysis, we find a variety of tempera- 
tures (from 0.2 to 0.7 keV) and Nh (from Galactic to 2 x 10 21 cm -2 ). Metal abundances for Ne, 
Mg, Si, and Fe vary dramatically throughout the ISM from sub-solar values (~0.2) up to several 
times solar. 

Subject headings: galaxies: peculiar — galaxies: individual(NGC4038/39) — galaxies: interac- 
tions — X-rays: galaxies — X-rays: ISM 



1. Introduction 

The Antennae, dynamically modeled by Toomre & Toomre (1972) and Barnes (1988), are the nearest 
pair of colliding galaxies involved in a major merger (D = 19 Mpc for H = 75 km s -1 Mpc -1 ). Hence, 
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this system provides a unique opportunity for getting the most detailed look possible at the consequences 
of a galaxy merger as evidenced by induced star formation and the conditions in the ISM. Each of the two 
colliding disks shows rings of giant H II regions and bright stellar knots with luminosities up to My <~ — 16 
(Rubin, Ford, & D'Odorico 1970), which are resolved with the Hubble Space Telescope (HST) into typically 
about a dozen young star clusters (Whitmore & Schweizer 1995; Whitmore et al. 1999). These knots 
coincide with the peaks of Ha, 2.2/i, and 6-cm radio-continuum emission (Amram et al. 1992; Stanford et 
al. 1990; Ncff & Ulvestad 2000), indicating an intensity of star formation in each single region exceeding 
that observed in 30 Doradus. CO aperture synthesis maps reveal major concentrations of molecular gas, 
including a ~2.7xl0 9 M Q concentration in the region where the two disks overlap (Stanford et al. 1990; 
Wilson et al. 2000). This concentration in the overlap region is coincident with the bright emission peaks 
seen in SCUBA 450^ and 850^ maps and ISOPHOT 60/x and 100^ observations (Haas et al. 2000), which 
show the simultaneous presence of warm dust (T « 30 K), usually observed in starbursts of luminous IR 
galaxies, and cold dust (T < 20 K) more typical for quiet galaxies and extremely dense clouds. A K-band 
study derives ages for the star clusters ranging from 4 to 13 Myr, and yields high values of extinction 
(Ay ~ 0.7 — 4.3 mag; Mengel et al. 2000, 2005). Many of the youngest star clusters lie in the CO-rich 
overlap region, where even higher extinction is present (Ay < 100 mag), implying that they must lie in front 
of most of the gas (Wilson et al. 2003). Some of them are found as far as 2 kpc from detectable CO regions, 
leading to the conclusion that either they formed from clouds less massive than 5 x 10 6 M Q or they have 
already destroyed their parent molecular cloud, or perhaps both. 

The presence of an abundant hot ISM in The Antennae was originally suggested by the first Einstein 
observations of this system (Fabbiano & Trinchieri 1983), and has since been confirmed with several major 
X-ray telescopes (ROSAT: Read, Ponman, & Wolstencroft 1995, Fabbiano, Schweizer, & Mackie 1997; and 
ASCA: Sansom et al. 1996). The ROSAT HRI image was used in a recent multiwavelength study of The 
Antennae by Zhang, Fall, & Whitmore (2001), which suggests feedback between star clusters and the inter- 
stellar medium (ISM). The first Chandra ACIS (Weisskopf, O' Dell & van Speybroeck 1996) observation of 
The Antennae in 1999 (Fabbiano, Zezas & Murray 2001; Fabbiano et al. 2003, hereafter F03) gave us for 
the first time a detailed look at this hot ISM, revealing a complex, diffuse and soft emission component re- 
sponsible for about half of the detected X-ray photons from the two merging galaxies. The spatial resolution 
of Chandra is at least 10 times superior to that of any other X-ray observatory, allowing us to resolve the 
emission on physical scales of ~75 pc for D=19 Mpc, and to detect and subtract individual point-like sources 
(most likely X-ray binaries; see Zezas et al. 2002). At the same time, ACIS allows us to study for the first 
time the X-ray spectral properties of these emission regions, providing additional important constraints on 
their nature. 

In the first detailed Chandra study of the hot ISM, F03, using the 1999 data, presented a "multi-color" 
X-ray image that suggests both extensive absorption by the dust in this system, peaking in the overlap region, 
as well as variations in the temperature of different emission regions of the hot ISM. These conclusions were 
confirmed by spectral fits to the data extracted from different regions, which suggested that at least two 
thermal components with temperatures of ~0.3 and ~0.7 keV are present in the luminous regions of the 
diffuse emission. The parameters derived for the hot ISM suggested that in the nuclear regions of the 
two colliding galaxies the mass of the hot ISM is ~10 5 -10 6 M Q , much smaller than that in molecular gas 
(Wilson et al. 2000), but that the thermal pressures of the two phases of the ISM are comparable, suggesting 
equilibrium. Comparisons with Ha, radio-countinuum (Neff & Ulvestad 2000), and HI data (Hibbard et al. 
2001) suggested a displacement of the diffuse X-ray emission on the northen side of the system, consistent 
with the effect of ram pressure from the cold ISM (H I) raining onto the hot gas in the northern star-forming 
arm. F03 did not make any conclusive statement on the metal abundance of the hot ISM. Metz et al. (2004) 
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analyzed the 1999 data together with the second and third observations, concentrating on regions of high 
star-formation rate, and found that super star clusters play a significant role in heating the ISM, but that 
the amount of hot interstellar gas is not directly proportional to the cluster mass. 

While this first Chandra data set demonstrated the richness of the ISM in The Antennae, the number 
of detected photons was such that the detailed small-scale morphology and spectral properties could not 
be fully explored. We now have completed a deep monitoring observing campaign of The Antennae with 
Chandra ACIS, that has produced a very detailed and rich data set. Our first look at these data revealed 
an intricate diffuse emission, with clear signatures of strong line emission in some regions, pointing to high 
metal abundances of the hot ISM (Fabbiano et al. 2004). In the present follow-up paper, we report a detailed 
analysis of this diffuse emission. The paper is organized as follows. In §2 we describe the Chandra ACIS- 
S observations of The Antennae and the procedures followed to screen and prepare the data for analysis. 
The creation of a broad-band mapped-color image of the diffuse emission in The Antennae is described 
in §3, while the preparation of a mapped-color line-strength image — useful for detecting differences in the 
spatial distribution of the metals — is described in §4. In §5 we present the spectral analysis performed after 
subdividing the distribution of diffuse emission into 21 regions, following the indications from the mapped- 
color images. We then discuss the adopted spectral models and the method used for computing the errors, 
present the results of the spectral fits, and compare our findings with previous results obtained by Sansom 
ct al. (1996) and F03. We summarize the results of our analysis in §6. In a companion paper (Baldi et al., 
submitted) we discuss the implications of these results for our understanding of the physical conditions and 
metal enrichment in the hot ISM. 



2. Observations and Data Preparation 

NGC 4038/39 were observed with Chandra ACIS-S seven times during the period between 1999 De- 
cember and 2002 November, for a total of ~411 ks (Table 1). The satellite telemetry was processed at the 
Chandra X-ray Center (CXC) with the Standard Data Processing (SDP) procedure, to correct for the motion 
of the satellite and to apply instrument calibrations. The data used in this work were processed in custom 
mode with the version R4CU5UPD6.5 of the SDP, to take advantage of improvements in processing software 
and calibration, in advance of data reprocessing. Verification of the data products showed no anomalies. 
In addition, we have reduced positional errors by using stars in the field detected in X-rays (see Zezas et 
al. 2002). This revised astrometry is consistent with absolute source positions with 1-er errors of 0'.'6. We 
have used the latest and best Chandra positional information for our comparison with other positional infor- 
mation. The data products were then analyzed with the CXC CIAO v3.0.1 software and XSPEC package. 
CIAO Data Model tools were used for data manipulation, such as screening out bad pixels and producing 
images in given energy bands. The latest release of the ACIS CCD calibration files (FEF) was used for the 
analysis. 

The data were screened to exclude the two "hot" columns present in ACIS-S3 at chip-x columns 512 
and 513 (at the b and c node boundaries). No screening to remove high background data was necessary: 
the light curve extracted from an area of 7 arcmin 2 around the observation aimpoint (excluding the point 
sources) showed a constant count rate of ^0.15 counts s _1 for every pointing but the last one (ObsID: 3041), 
during which the background rate rose to ~ 0.45 counts s _1 for about 25 ks. However, this effect is relatively 
mild (at worst three times the "quiet" background) and on balance the removal of 1/3 of this observation 
would be more detrimental to the data analysis than the slightly increased uncertainties deriving from a 
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larger background level. 

We applied corrections to account for the temporal evolution of the ACIS spectral response, caused 
by changes in the Charge Transfer Inefficiency (CTI) and in the electronic gain of the CCDs. This is a 
sizeable effect in our case, which we corrected for by applying the time-dependent gain correction described 
in http:/ /cxc. harvard.edu/contrib/alexey/ tgain/tgain.html. Figure 1 shows a spectrum of a diffuse emission 
region, plotted before and after the gain correction. A shifting in the emission lines present in the spectrum 
is evident, with consequences on the quality of the spectral fits (see residual plots in Fig. 1). 
The entire data set was then co-added to produce the image of the central part of The Antennae shown in 
Figure 2. 

3. The mapped-color image of the diffuse ISM 

Figure 2 shows that complexly structured diffuse X-ray emission is clearly present in The Antennae. 
Previous work on the 1999 observation alone has shown that this diffuse emission can be related to a hot 
ISM with varying spatial and spectral properties (Fabbiano et al. 2001; F03). However, our much deeper 
data show a richness of detail in this diffuse emission that was not visible in the 1999 data alone. 

To study this diffuse emission, we followed the procedure of F03 to first subtract the point-like X-ray 
sources from the data, fill the resulting holes by interpolating the surrounding diffuse emission, and then 
produce a mapped-color image of the hot ISM. First, from the data we extracted images in three different 
energy bands: 0.3-0.65 keV, 0.65-1.5 keV and 1.5-6 keV. The middle energy band was chosen to encompass 
the emission from the Fe-L blend. The 0.3 keV low-energy boundary is the lowest energy within which ACIS 
is well calibrated. The upper energy boundary was set at 6 keV because at energies >6 keV the diffuse 
emission becomes background-dominated. Then, we used the CIAO tool wavdetect to generate a list of 
sources detected at a threshold of 3a in the co-added exposure, using data in the total 0.3-6 keV energy 
band within the central 2' of the image. From this source list, we removed the extended sources listed in 
Zezas et al. (2002) plus a few additional slightly-extended sources identified by visual inspection of the data. 
This culling left a total of 127 point-like sources, the positions of which are shown in Figure 2. 

We then used the CIAO tool dmfilth to subtract these 127 sources from the data. This tool cuts out 
the point sources and fills in the holes corresponding to the source positions with values interpolated from 
the surrounding background regions. We adopted as source regions the 3a best fit ellipses computed by 
wavdetect. As background regions we adopted elliptical annuli having the source region ellipse as inner 
boundary and an ellipse with the same orientation, but with double size semi-minor and semi-major axes, 
as outer boundary. Areas in the background regions where unrelated sources happened to fall were excluded 
from the calculation of the background. We used an interpolation technique in which pixel values in the 
source regions are sampled from the distribution of pixel values in the background regions (dist option in 
dmfilth), since the Poisson interpolation in the current version of dmfilth does not work for low background 
regions (< 0.10 cts/pixel) such as those in the outer parts of The Antennae. 

The three diffuse-emission images at different energies were then adaptively smoothed, using csmooth. 
This process preserves the high-significance small-scale details of the image, while enhancing more extended 
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features of low surface brightness. The same smoothing scales were used for the three images, based on those 
computed during the smoothing of the 0.65-1.5 keV image, which is the one with the highest signal-to-noise 
ratio. For the smoothing, we used a Gaussian kernel of varying size (from 1 to 40 pixels). The maximum 
significance was set to 5a, while the maximum smoothing scale corresponds to a minimum significance of 
~2.4<j for the 0.6-1.5 keV image. The resulting smoothed images were then corrected for exposure (i.e., 
flat-fielded) by dividing them by the appropriate exposure maps, smoothed using the same scale as the 
image. The exposure maps were generated for a thermal emission model (Raymond & Smith 1977) with a 
temperature of 0.8 keV. The resulting exposure-corrected and adaptively-smoothed images, shown in Fig- 
ure 3, were then combined (using the tool dmimg2jpg) to create the "mapped-color" image of the diffuse 
emission, shown in Figure 4. 



In Figure 4, red denotes softer X-ray emission regions (0.3-0.65 keV), while blue denotes harder emission 
regions (1.5-6 keV). Green identifies emission in the middle band (0.65-1.5 keV). In the following discussion, 
we will refer to the counts in these bands as S (for soft) , M (medium) , and H (hard) . The color bars in this 
figure quantify the distribution of the pixel values (from minimum to maximum) in each band. 



By comparison with the mapped-color image obtained from the first dataset (ObsID: 0315; F03), the 
new, deeper image shows an impressive level of detail. Spatial and spectral structures are visible on different 
spatial scales, from the smaller-scale clumps in the inner parts of the optical merging galaxies to the larger- 
scale regions of low surface brightness . 



Since we have three color-coded energy bands, the resulting pixel color will be ideally distributed in a 
tridimensional "color cube" whose axis values are proportional to the number of counts in a given band (see 
right panel of Figure 5). We have "sliced" this cube to obtain a two-dimensional representation that relates 
the pixel color to the hardness ratios HRl and HR2, defined as: 

HR1 = ^; HR2=^-. 
M+S H+M 

Each of these slices represents a plane with the same value of hard-band counts (see Figure 5), with soft-band 
counts increasing along the horizontal axis and medium-band counts increasing vertically. For this chosen 
value of hard counts, each color in the slice can be uniquely related to a value of HRl. These values are 
written in the bottom slice. The same values will apply to each slice. HR2 does not depend on the amount 
of soft counts, so we display a column of HR2 values for each slice, along the "medium-band" axis. 



It is possible to also give a physical meaning to these numbers, computing the values of HRl and HR2 
in two cases: an absorbed Raymond & Smith (1977) thermal model and an absorbed power-law model. For 
the thermal model we vary the temperature from 0.2 keV to 0.8 keV and the metallicity from 0.1 solar to 2 
times solar. In the power-law case the photon index T ranges from 1 to 3 and the intrinsic absorption from 
10 20 to 10 21 " 5 cm~ 2 . In both cases we considered also a Galactic Nh component fixed at the line-of-sight 
value toward The Antennae (~ 3.4 x 10 20 cm~ 2 ; Stark et al. 1992). The values of HRl and HR2 are plotted 
in the X-ray color-color diagram shown in Figure 6, and they can be directly related to the labels in Figure 5. 
The red grid refers to the thermal model, while the blue grid represents the power-law model. It is likely 
that a simple single-temperature thermal model or an absorbed power-law model may not be a good fit to 
the data. However, the figure gives an indication of the characteristics of the X-ray emission from different 
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regions. 

For example, the blue region on the bottom-left of Fig. 4 corresponds roughly to HR1 ~ 0.7 and 
HR2 ~ —0.5. Looking at Figure 6, these hardness ratios — marked by a black square — are located at a point 
in the diagram between the absorbed-power-law model and the purely thermal model, possibly indicating 
the presence of unresolved binary emission in this region. On the other hand the orange region on the right 
of Fig. 4 corresponds roughly to HR1 ~ 0.5 and HR2 <~ —0.9 (black triangle in Figure 6). These values are 
typical of purely thermal emission. 



4. The Line-Strength Map 

As Figure 1 illustrates, the spectra extracted from different regions of the diffuse emission frequently 
show emission lines. In particular, O+Fe+Ne emission (0.6-1.16 keV) is often present, together with Mg-XI 
(1.27-1.38 keV) and Si-XIII (1.75-1.95 keV) lines. Using these features, we generated another color map 
(see Fabbiano et al. 2004) that gives a qualitative representation of the presence of emission lines from Iron, 
Silicon, and Magnesium in different regions of the hot ISM. This map has proven very useful for identifying 
regions of strong apparent line emission and for selecting spectrally similar regions for a proper spectral 
analysis. 

The line-strength map, shown in Fig. 7a, was constructed as follows. First, we created images in three 
bands encompassing emission from O+Fe+Ne (0.6-1.16 keV), the Mg-XI line (1.27-1.38 keV), and the Si- 
XIII line (1.75-1.95 keV). We also created an image including emission in the 1.4-1.65 keV and 2.05-3.05 keV 
bands, which was used for continuum subtraction. The line and the continuum bands were determined from 
spectra of the diffuse emission in order to have an optimum line/continuum ratio and minimal contamination 
by emission lines, respectively. In each band we excluded all point-like sources and interpolated over the 
holes as described in § 3. The O+Fe+Ne band image, which contains a much larger number of counts than 
the other line-band images, was used to determine the smoothing scales, and was adaptively smoothed to 
a significance between 3a and 5a. The same smoothing scales were then also applied to the other line and 
continuum images. In order to subtract the appropriate continuum from each line image, the continuum 
image was rescaled to represent the band continuum in the given line. To calculate the weighting factors, 
we extracted the spectrum of the entire diffuse emission of The Antennae and fitted the continuum (after 
excluding the energy ranges with line emission) with a two-component Bremsstrahlung model. Based on 
this model we estimated the relative intensity of the continuum in our continuum band and the three line 
bands. Each line image was continuum subtracted after weighting the continuum image by these weighting 
factors, and the resulting images were combined to create a mapped-color line-strength image (Figure 7a). 
In this image, red represents line emission from O+Fe+Ne, green from Magnesium, and blue from Silicon, 
respectively. 

Figure 7a shows that the hot ISM of The Antennae has a complex emission-line structure. Silicon 
emission (blue) appears pervasive in most of the hot ISM, with Magnesium lines (green) also present in 
various relative amounts and O+Fe+Ne emission (red) prevalent in some regions. Other regions (black) 
show only very weak or no line emission. 
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5. Spectral Analysis 

The mapped-color X-ray image (Figure 4) and the line-strength map (Figure 7a) demonstrate that the 
diffuse emission of the hot ISM is rich in spatial and spectral features, and that the line emission in the X-ray 
spectra is pervasive and varied. These figures suggest a temperature structure in the hot ISM, possibly the 
effect of varying absorption on the emission (see also F03), and of its metal enrichment. 

5.1. Region Selection and Extraction of the X-ray Spectra 

To put these results on a more rigorous and quantitative footing, we have performed a spectral analysis 
of different regions of the hot ISM, using Figures 4 and 7a as guides for the selection of these regions. We first 
identified 17 separate regions for analysis, based on the morphology of the diffuse X-ray emission (Fig. 7a). 
However, the line-strength map suggested marked abundance differences within some of these regions. In 
particular, Figure 7a shows that regions 4, 6, 8, and 12 have a complex metallicity structure. For example, 
in Region 4 we can identify a metal-rich and a metal-poor area; in regions 6, 8, and 12 we see clearly defined 
sub-areas with yellow (O+Fe+Ne and Mg) and red (O+Fe+Ne) emission. These four regions also have a very 
high number of detected counts, making it possible for us to subdivide them further on the basis of their line 
features. We subdivided each of these four regions into two parts, as shown in Figure 7b. Although Figure 7a 
suggests spatially complex metal distributions also in other regions (5, 13, 14 and 15), these regions have 
too few counts to be subdivided further. The described process identified a total of 21 regions for separate 
spectral analysis. 

From these 21 regions, we excluded all the point source 3a areas described in § 3 and shown in Figure 2. 
We then extracted background counts representative of the field from three circular source-free regions well 
outside the galaxies, but within the same CCD (ACIS-S3) in all seven observations. These background 
regions covered a total area of ^4 arcmin 2 . 

All of the data sets under analysis, with the exclusion of the first and fifth (ObsID: 315 and 3044, 
respectively), were obtained with the same ACIS-S configuration and could then be treated as a single 
observation. We therefore merged these data into a single event file by using the CIAO script merge-all. 
ObsID 315 and 3044, however, needed individual handling. ObsID 315 was obtained in 1999 December, when 
the low-energy response of ACIS-S was significantly higher (A e f / ~ 500 cm 2 , QE <~ 65% at 1 keV) than at 
the time of the other six observations (2001 December to 2002 November). ObsID 3044 was obtained with 
a different ACIS-S configuration (different SIM-Z). 

We extracted spectra separately from the merged 5-ObsID file and from 315 and 3044, using the diffuse 
emission and background areas described above. In all cases, the CIAO script acisspec was used to extract 
the region spectra. This script invokes the tool dmextract to extract the spectrum and then uses the tools 
mkwrmf and mkwarf to create area-weighted Response Matrix Files (RMF) and Ancillary Response Files 
(ARF) for each region, since they span more than one FITS Embedded Function (FEF) region. The spectra 
derived from the three different event files, both from source and background regions, were then summed 
using the FTOOL mathpha. The response files were combined with their appropriate weights, using the 
FTOOLS addrmf and addarf. The resulting spectra were then grouped so as to have at least 20 counts 
per energy bin before background subtraction. 

For each of the 21 regions used in the spectral analysis the net counts (background subtracted) in three 
different energy bands and in the total 0.3-6.0 keV band are listed in Table 2. 
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5.2. Spectral Fit Method 

We analyzed the extracted spectra with XSPEC vll.2.0 (Arnaud 1996), using combinations of optically- 
thin thermal emission and power-law components to fit the data. We used the Astrophysical Plasma Emission 
Code (APEC) thermal-emission model (Smith ct al. 2001) to represent the thermal emission. This model 
is based on the original Raymond & Smith (f977) code, but has been totally revamped to exploit modern 
computing capabilities and the wealth of accurate atomic data. 

The real nature of a soft X-ray emitting gas, recently heated as in star-forming galaxies, is clearly 
that of a multi-temperature plasma with very complex differential emission- measure distributions (see, e.g., 
Strickland & Stevens 2000). Since the gas is continuously being heated by SN and star winds, it may also 
not be in ionization equilibrium. However, modeling our low S/N spectra with complex models having more 
than two temperatures or not assuming ionization equilibrium would require too many free parameters to 
obtain any meaningful information from them. For example, although fitting higher S/N spectra of diffuse 
X-ray emission in nearby starburst galaxies with two-temperature models (using combinations of several 
equilibrium and non-equilibrium models) frequently leads to better fits than fitting with a single-temperature 
model (e.g., NGC 253: Strickland ct al. 2002), they almost always lead to unrcalistically low abundances, as 
previously noted by Metz et al.(2004). Therefore, we considered also single-temperature models to fit our 
data, together with two-temperature models, to allow us to compare the results obtained from both. The 
four models that we used, listed in order of increasing complexity, are: 

• absorbed one-temperature APEC model (XSPEC model: wabs(wabs(vapec))); 

• absorbed one-temperature APEC model plus a power-law component (XSPEC model: wabs(wabs(vapec+pow))); 

• absorbed two-temperature APEC model (XSPEC model: wabs(wabs(vapec+vapec)J); and 

• absorbed two-temperature APEC model plus a power-law component (XSPEC model: wabs(wabs(vapec+vapec+pow))). 

Two different absorption components were used in all cases: one fixed at the Galactic line-of-sight value 
toward The Antennae {N H = 3.4 x 10 20 cm" 2 ; Stark et al. 1992) and the other one left free to vary in order 
to model the intrinsic absorption within the two galaxies. In the two models where a power-law component 
is included we tried also to fit the Nh absorption independently for the power-law and thermal components. 
Although the Nh of the power-law component resulted in a value different from the Nh of the thermal 
component(s) (on average ~ 10 21 cm~ 2 instead of ~ 10 20 cm -2 ), the best fit parameters of the thermal 
component (mainly kT and Z) did not vary significantly and were fully consistent within the errors. The 
abundances of Neon, Magnesium, Silicon, and Iron were all left free to vary. The other elements (e.g., 
Carbon, Nitrogen, and Sulfur) arc poorly constrained because of gain uncertainties and sensitivity to Nh- 
Thus we adopted solar abundances for these elements. We tried to thaw the Oxygen abundance, obtaining 
no clear improvements in the fits and no significant variations in the best-fit parameters of kT and of Z in 
the other elements. Moreover, the poorly costrained Oxygen abundance values were always consistent with 
solar values. Therefore, we decided to adopt a solar abundance also for this element. 

We introduced a power-law component into our models to allow for the possible residual presence of un- 
resolved X-ray binaries too faint to be detected individually. Indeed, from an analysis of the entire emission 
throughout The Antennae we found that the point sources contribute up to ~85% of the emission in the 2-5 
keV range and are essentially the only contributors above 5 keV. The photon index value was at first frozen 
at T = 1.88, a value obtained by fitting the 2-10 keV spectrum of the entire emission from The Antennae, 
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including point sources. Afterwards F was left free to vary to fit the data in order to test for the presence of 
substantial discrepancies between the average power law determined from the galaxies' global emission and 
possible power-law components in individual regions. Wc compared the results of the two fits (frozen and 
thawed T) for the same region using an F-test, and considered the thawed T model to be a better description 
of the data if p, the probability that the improvement in the x 2 statistics is due by chance, was less than 5%. 

Fhc resulting best-fit temperature(s) kT, absorbing column density Nh, and power-law photon index 
r (if a power-law component is necessary) for the 21 regions studied are summarized in Table 3, together 
with the abundances of Neon, Magnesium, Silicon, and Iron relative to the solar values. The last column of 
Table 4 gives the "best-fit model". The errors are quoted at la for one interesting parameters. The values 
listed in the table are obtained from the "best-fit model", which is defined as the simpler model (i.e. with 
a smaller number of parameters to be fitted) with x 2 e d < T4. This value of x 2 e d represents a fair boundary 
between an acceptable fit and a bad fit for low S/N spectra (0.3-6 keV counts in the range between ~ 750 
and ~ 6500). This means that we do not adopt a two-tcmpcraturc model if an acceptable fit is already 
obtained with only one temperature. We understand that a singlc-tcmperature and even a two-temperature 
model are approximations to the real physical nature of the X-ray emitting gas, which is likely to be more 
complex. However, given the limitations of our data we followed the customary data analysis practice of 
choosing the simpler model of two if both models give acceptable representations of the data. Indeed, we 
found that only one region does not yield an acceptable fit using a single-temperature spectral model (Region 
4b). Possible systematic biases depending on the model choice are taken into account as explained in § 5.2.1 
below. The only region where none of our models yielded a Xred < T4 (although with only 751 counts in 
the 0.3-6 keV band) is Region 13, for which we chose the model with the lowest xled- The spectra of the 
21 regions, together with the best-fit models and the data-to-best-fit-model ratios, are shown in Figures 8, 
9, 10, and 11. 



5.2.1. Error computation 

Although quoting the errors at la for one interesting parameter (as the errors quoted in Table 3) is a 
common practice in X-ray spectral analysis, it may lead to an underestimation of the true uncertainties not 
only because of possible systematics but also in the case that some of the fitting parameters are correlated. 
As previously noted by Martin, Kobulnicky & Heckman (2002) a degeneracy between metallicity and tem- 
perature may exist in the spectral models, and it is very difficult to solve such ambiguity from low resolution 
CCD spectra like the ones we are dealing with. Indeed we found such a degeneracy at least in some of our 
regions. Although it will not be possible to solve completely this problem from the spectral analysis of ACIS 
spectra alone, it is still possible to adopt a very conservative approach to take into account both systematics 
and correlations between fitting parameters. This approach mainly consists in calculating the errors for the 
best-fit parameters at la for all interesting parameters. The errors on kT, Nh, and T were calculated from 
the x 2 statistics considering all free parameters in the "best-fit model" as interesting. Since the spectral 
regions occupied by Ne, Mg, and Si lines are quite distinct, the errors on the abundances for these elements 
were computed using a different method to reduce the number of interesting parameters. We first attempted 
to link together all the a elements using the type II SNe ratio from the compilation of average stellar yields 
from Nagataki & Sato (1998), in order to reduce the number of free parameters by two. However, we were 
unable to obtain acceptable fits for any of the regions, suggesting that the elemental ratios are not entirely 
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consistent with type II SN yield models. Thus we tried an alternative approach for each individual a-element 
line, ignoring from the fit the part of the spectrum relative to the other two lines (e.g. for Mg we ignored 
the part of the spectrum occupied by the lines of Ne and Si) and freezing their abundances to their best-fit 
values. Therefore, in the error computation of the a-elements we have two free parameters less than in the 
computation of kT, Nh, and T. For deriving the Fe abundance we ignored the part of the spectrum relative 
to all three a elements (freezing them to their best- fit values) , obtaining three less free parameters than the 
original model. The number of interesting parameters was set equal to the number of free parameters also 
for the abundance error estimation. Note that this is a very conservative approach for the determination 
of statistical errors, since the fitting parameters are not completely independent from each other, so that 
the number of truly interesting parameters is likely to be lower than the one adopted. However, statistical 
errors are not the only uncertainties present in our measurements, and systematic biases may arise also 
from the spectral model choice, as noted before in § 5.2. Although our method cannot estimate exactly 
any possible systematics, it is conservative enough that we are confident about the parameter confidence 
ranges so derived. The best-fit parameters and the errors computed using this method are shown in Table 4. 
Comparing this table with Table 3, it is clear that the abundances measured using the method described 
above are consistent with the measures performed fitting the whole spectrum. The only noticeable difference 
comes from the la confidence ranges which, of course, are larger in Table 4. In the rest of the paper we will 
use the confidence ranges of Table 4. 

The x 2 parameter space is often asymmetric. Therefore, we performed an analysis of the behaviour of 
X 2 with each free parameter (using the steppar command in XSPEC), to explore the presence of possible 
secondary x 2 minima in the confidence range estimation. Figures 12-16 show some examples of the most 
complicated features observed in the x 2 parameter space when varying the temperature (kT) and the abun- 
dances of Ne, Mg, Si and Fe, respectively. As can be seen in Figure 12, one of the temperatures of Region 
4b is unconstrained even at the la level. Three other regions (8b, 14 and 17) show two minima in x 2 space. 
For them we adopted la confidence ranges encompassing both minima even when the Ax 2 value between 
them exceeded the la confidence level, as happens in Regions 8b and 14. For abundances the presence of 
secondary minima is more rare: we found a secondary minimum only in Region 8b (Si and Fe; Figs. 15 and 
16) and in Region 4b (Fe; Fig. 16). 

5.2.2. Emission-line equivalent width 

While spectral fits are the standard way to obtain elemental abundances from X-ray data, we also wish 
to obtain abundance limits that are relatively model-independent and that allow us to evaluate the possible 
effects of departures from ionization equilibrium (NEI) . Since in our spectra the lines of Mg and Si are not 
confused with other elements, they allow us to directly measure their equivalent widths (hereafter EW). 
For each region we considered the spectrum in the whole 0.3-6 keV range and used the best-fit model (as 
defined in § 5.2) for the continuum. We artificially set to zero the abundance of the element in consideration, 
adding a gaussian line to the model. Both the normalization and the energy of the gaussian line were left 
free to vary in the fit. From the best fit so obtained we measured the EW of the line using the XSPEC task 
eqwidth, using only the thermal component (s) as continuum emission. With eqwidth we performed a Monte 
Carlo sampling of the parameter values distribution, estimating the errors on the EW (at la) independently 
from the x 2 statistics. The relation between the measured EW and the measured abundance for Mg and Si 
are shown in Figure 17. We find a clear correlation between the two parameters for Mg; a weaker correlation 
is suggested by the Si plot, however we must remember that the continuum is much less well defined for this 
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line, increasing the uncertainties. 

The metal abundances measured from the spectral fits can be directly related to the EWs by the 
expression: 



where Z is the abundance of the element and a and a are two parameters essentially dependent only on 
the temperature of a hot X-ray emitting gas in a static ionization equilibrium. The denominator comes 
about because the continuum is produced by bremsstrahlung and the recombination of lower atomic number 
elements, while the numerator is simply proportional to the abundance of the element in question. The 
temperature dependence arises mostly from the ionization fractions of the H-like and He-like ions that 
produce the emission lines. The term aL in the denominator accounts for the contributions of all the 
elements other than H and He to the continuum through bremsstrahlung, recombination, and two-photon 
emission assuming that all the elements scale together. Thus at very high Z the EW will reach an asymptotic 
limit, which can only be exceeded if the clement in question is preferentially enhanced relative to the elements 
that contribute to the nearby continuum or if NEI increases the equivalent width. In order to assess the likely 
effects of NEI we computed equivalent widths from shock wave models using an updated version of the code 
described in Raymond (1979) and Cox and Raymond (1985). While the EW can be far from the equilibrium 
values just behind the shock, it tends to be dominated by emission from the region that is not too far from 
equilibrium. We computed models for shocks in the range 400-900 km/s, corresponding to a temperature 
immediately behind the shock between 0.19 and 0.95 keV. Note that in the latter case the temperature that 
would be derived from the X-ray continuum would be somewhat lower because of contributions from gas in 
the cooling region. 

Using the values of a and a computed for the two different cases examined (thermal plasma and shock 
waves) we derived the expected EW corresponding to an abundance of jZq, Z q and 5Z Q for both Mg and 
Si. Figure 18 shows the relation between the EW of Mg and the EW of Si for all the regions where we 
were able to obtain a meaningful measure of both. The regions of the diagram relative to abundances of 
\Zq, Zq and 5Zq, as a function of the temperature are also plotted. If we consider the simple equilibrium 
thermal model (Fig. 18, left) the majority of the points are located in the region corresponding to solar 
abundance, but there are regions located in the upper-right portion of the diagram as well. Region 5 has 
EW not consistent with Z < 5Z@ for both elements, in agreement with the results of the spectral analysis 
where this region shows the highest abundances of Mg and Si among all the spectra analyzed. Also Regions 
7, 8b, 15, and 17 show high values of the EW, but with large errors, consistently with what we observe in 
the abundances measured from the spectral analysis of these regions. The situation is even more dramatic if 
we consider the shock model (Fig. 18, right), where the values of Mg for the regions listed above are clearly 
not consistent with Z < 5Z Q (as the spectral analysis would suggest), and the Silicon EW for Region 7 is 
also not consistent with Z < 5Z Q . Notice that the values of the EW we derived from the abundances at 
the different temperatures were computed assuming that all the elements scale together at the same time. 
Thus the continuum due to lower Z elements contributes at the energies of the emission lines. If Mg or Si 
are enhanced relative to O in particular, the equivalent width expected could be higher. That would require 
that the enrichment in The Antennae is done by supernovae with a different mass distribution than those 
that enriched the solar neighborhood (e.g. Woosley & Weaver 1995). This is indeed likely for very young 
star forming regions, in which only the most massive stars will have completed their evolution. We note that 
departures from ionization equilibrium other than those of the shock models are possible if the plasma is 
rapidly heated or cooled. In general they do not raise the EW very much above the equilibrium values except 
in a very strongly overionized, cool plasma (e.g. Breitschwerdt & Schmutzler 1994). While it is possible to 
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produce such a plasma by rapid adiabatic expansion, the low emissivity makes it difficult for such a plasma 
to make a significant contribution to the spectrum. 



5.3. Comparison with results from the first Chandra observation 

Almost all the regions (20 out of 21) did not require the presence of a second thermal component in the 
model to fit the data, and half of them (11 out of 21) did not need the presence of a power-law component to 
fit the data. In the cases where the power- law component was needed, the fixed T = 1.88 models were found 
to yield a better fit than a thawed T model in eight regions out of 10. Only Region 4b required two thermal 
components with different kT. In this region, we could not constrain the low kT value (although we found a 
minimum in the x 2 statistics at kT = 0.20) and we could place only a lower limit (kT > 0.54) on the higher 
temperature component. Almost all the single-temperature regions required a temperature of kT ~ 0.6 keV 
to fit the data, with regions 5, 8b, and 14 (kT ~ 0.3 keV) the only exceptions to this rule. Considering all 
the regions analyzed, the temperatures (both in the single- and double-tcmpcrature cases) range from ^0.2 
to ~ 0.8 keV (considering also the errors), a range fully consistent with the values found in F03. 

F03, however, could not reach any conclusions about the metal abundances of the ISM because of the 
shorter integration time of their data set (only ObsID 315). In contrast, our analysis yields metal abundances 
that have acceptable constraints in the majority of the regions. This represents the most striking new result 
of our spectral analysis: the abundances of Silicon and Magnesium cover a wide range of values from the 
very low values found in regions like Region 2, 6a, 12a, and 12b (Z ~ 0.2Z Q ) to the huge values observed in 
Region 5 (Z > 20Z Q ). 

The intrinsic absorption throughout the hot ISM is generally low, with typical Nh ~ 10 20 cm~ 2 and 
often consistent with zero. The exceptions are the southern nucleus (regions 8a and 8b) and Region 7 (the 
dark-blue area in Figure 4) , which is significantly obscured also in the optical. Region 7 corresponds to the 
Overlap Region, where the most active star formation is now occurring (Mirabel et al. 1998; Wilson et al. 
2000; Zhang et al. 2001). These regions show intrinsic absorption corresponding to Nh ~ 1 — 2 x 10 21 cm~ 2 . 

A direct comparison with the results of the spectral analysis performed on ObsID 315 in F03 can be 
done for the four emission regions identified and analyzed individually in that paper: the northern nucleus 
(roughly corresponding to our regions 14 and 15), the southern nucleus (regions 8a, 8b, and 9), the hot spot 
Rl (regions 12a and 12b), and the hot spot R2 (Region 13). 

In the northern nucleus F03 found kT ~ 0.6 keV, which is fully consistent with the temperatures we 
found in regions 14 and 15 (0.37±g; | and 0.59±°;^ keV, respectively) . The best-fit abundances given in F03 
(~0.15 solar) are lower than our solar or slightly sub-solar abundances observed in Region 15 (in Region 14 
abundances are unconstrained); however, the large errors encompass our results. 

The southern nucleus in F03 is well fitted by a two-temperature model, with temperatures (0.35 and 
0.82 keV) consistent within the errors with the two peaks in % 2 -space observed for the temperatures of our 
Region 8b (0.33 and 0.60 keV). However, in regions 8a and 9 we find a single temperature gas at ~0.6 keV. 
Also in this case the abundances in F03 are very low (~0.2 solar), in contrast to abundances consistent with 
solar in all the three regions we examined. In this case as well, though, the abundances determined by F03 
are poorly constrained and consistent with our results. 
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In F03 the hot-spot Rl was fitted with a single-temperature model at ~0.4 keV, while the temperatures 
we find are somewhat higher, with 0.601q'q9 keV and 0.61 ± 0.06 keV for regions 12a and 12b, respectively. 
F03 determined a value of ^0.2 solar for the abundance, which is fully consistent with the low values we find 
in both Region 12a and Region 12b. 

For the hot spot R2 the only comparison possible is for the temperature (F03 could not constrain the 
abundance), which was ^0.4 keV in F03, somewhat lower than the value in Region 13 (£:T = 0.58±°;°f keV). 
However, the large errors in our new temperature determination (due to a secondary minimum in the x 2 , as 
described in § 5.2.1) make these two results consistent. 

The lower metal-abundance values suggested by the F03 analysis could be the result of blending of 
complex emission regions. Note, however, that in F03 all elements were kept locked at the solar ratio, while 
we have here fitted the abundances of Fe, Mg, Ne, and Si individually. 

Some of the 21 regions we have chosen are similar to, though not quite the same as, several of those 
analyzed by Metz et al. (2004). Our regions 2, 4, 5, 6, 12, 13 14 and 15 roughly correspond to their regions 
D03, D06, D07, D08, D01, D02, D04 and D05, respectively. In general, our and their spectral results are 
similar, except that the higher quality of our spectra permits at least a coarse measurement of elemental 
abundances in contrast with their assumption of solar abundances throughout the hot ISM. 

5.4. Comparison with ASCA Results 

The high metal abundances suggested by our spectral analysis in some regions of The Antennae are at 
odds with a previous claim, based on ASCA CCD spectra, of an overall extremely low abundance of heavy 
elements in this system (^0.1 the solar value, Sansom et al. 1996). Since the spectral resolution of ASCA 
was comparable to that of Chandra ACIS, we explore the possibility that — besides the lower signal-to-noisc 
ratio of the previous data — poor spatial resolution may have played a role in the low abundances reported 
from the analysis of the ASCA observations. 

Given the large (~ 2 arcmin) ASCA beam, the entire emission from The Antennae was used to derive a 
spectrum. To make a fair comparison with the ASCA data, we extracted from our Chandra data a spectrum 
of all 21 regions together, leaving the point sources (which could not be excluded in the ASCA observation) 
in the extraction region. We used the procedure described in § 5.1 to extract the spectra and compute 
response matrices individually for the first (December 1999) and the fifth (July 2001) observation, and from 
a merged event file for the other five observations combined. The background was taken from the same areas 
used in the case of the individual regions' analysis. The combined spectra and response matrices were used in 
XSPEC. The spectrum was fitted using a model with two additive components: a Raymond-Smith thermal 
model (for the diffuse emission) and a power law (for the point sources). To repeat the ASCA analysis we 
assumed a solar abundance ratio between the various elements, and fitted the metal abundance with just one 
parameter. We considered two different components for the absorption, one fixed at the Galactic value and 
the other left free to vary to measure the intrinsic absorption within The Antennae. The quality of the fit we 
obtained is very poor (xl ~ 2.8), not surprisingly given the complexity of the temperature and metallicity 
distributions throughout The Antennae. However, the best-fit model has a temperature of kT = 0.59 ± 0.03 
keV and a metallicity of O.lSt^Qg times the solar value, fully consistent with the ASCA results (Sansom 
et al. 1996). This demonstrates that the analysis of separate "clean" regions, possible because of the sub- 
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arcsecond resolution of Chandra, is needed to detect the metal lines in a complex hot ISM with CCD spectral 
resolution. 

The above conclusion is reinforced by the intercomparison of the spectral-fit results for different subsec- 
tions of morphologically selected regions (e.g., regions 5 and 12b), where the line-strength map of Figure 7a 
suggested different metal contents in the spectra. The results of these fits are compiled in Table 5 and clearly 
present intermediate metal abundances, typically lower than the peak value. We note that while the spatial 
resolution of Chandra surpasses that of any other X-ray telescope it is still finite, and while our exposure of 
The Antennae is deep the collecting area of Chandra is limited. Therefore, our results on metal abundances 
by necessity still suffer to some degree from an "averaging" effect over spectrally complex regions. As such, 
they may still represent lower limits on the local metal enrichment in some regions. 

This conclusion also agrees with unrelated recent results that demonstrate how fitting thermal plasmas 
with simplified (single-temperature) emission models may underestimate metal abundances (e.g., Molendi & 
Gastaldcllo 2002; Kim & Fabbiano 2004; see Fabbiano 1995 and Buote & Fabian 1998 for earlier discussions 
of this effect). 

6. Summary and Conclusions 

We have performed a detailed study of the X-ray properties of the diffuse emission of The Antennae 
(NGC 4038/39), analyzing the entire 411 ks exposure obtained as part of our monitoring observing campaign 
of this system with Chandra ACIS. Confirming the results of F03, which were based on the first of the seven 
observations used in the present study, we report a spatially and spectrally complex hot ISM. Thanks to our 
deep data, we also detect clear, spatially variable emission lines of Ne, Mg, and Si, in addition to the Fe-L 
blend, which indicate super-solar metal abundances in a few regions of the hot ISM. 

In summary: 

(1) We derive an X-ray broad-band color image of the hot ISM, as well as a color map displaying the 
apparent strength of line emission from O+Fe+Ne, Mg-XI, and Si-XIII. Both maps demonstrate complexity 
and have been used to guide our selection of 21 separate spectral extraction regions for detailed analysis. 

(2) Spectral models combining optically thin thermal-plasma and power-law emission can be fitted to 
our data with the presence, in almost all regions, of a single thermal component with temperatures ranging 
from 0.2 to 0.7 keV. Significant variations of Ah throughout the ISM are detected: column densities range 
from values consistent with the Galactic foreground absorption to higher values of ^2 x 10 21 cm~ 2 , as in the 
southern nucleus and in Region 7, the actively star-forming Overlap Region (Mirabel et al. 1998; Wilson et 
al. 2000; Zhang et al. 2001). A power-law component is present in half of the regions, showing values of T 
consistent with emission from unresolved X-ray binaries. The study of this power-law emission will be the 
subject of a forthcoming paper. 

(3) Fitting the Fe-L, Ne-IX, Mg-XI and Si-XIII emission, we find significant metal enrichment of the 
hot ISM. Metal abundances are generally consistent with solar, but reach extremes of subsolar and ^20-30 
solar in a few regions. In supporting evidence, the EW of the emission lines of Mg and Si correlates well 
with the abundances measured. Previous reports (from ASCA) of subsolar abundances in the hot ISM of 
the Antennae can be explained with an " averaging" effect due to the inclusion of the entire emission (diffuse 
from different regions and point sources) in a single spectrum. 
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The implication of these results for the physical state of the hot ISM and SN enrichment scenarios are 
discussed in a companion paper (Baldi et al., submitted). 
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Fig. 1. — Chandra ACIS-S spectrum of a prominent emission line region of the Antennae hot ISM, before 
applying the tgain correction (top) and after applying the tgain correction (bottom). The spectrum and 
the ratio with respect to a thermal+power-law model (XSPEC: wabs(vapec+powerlaw)) are plotted. 
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Fig. 2. — 0.3-6 keV image of the central part of The Antennae. The 3cr ellipses of the point sources detected 
by CIAO wavdetect (and reported in Zezas et al. 2002) are plotted in green. 



Fig. 3. — Adaptively smoothed images of The Antennae after source removal in the 0.3-0.65 keV band (red), 
0.65-1.5 keV band (green) and 1.5-6 keV band (blue). 



Fig. 4. — Mapped-color smoothed image of X-ray emission in The Antennae. The red color channel indicates 
emission from the 0.3-0.65 keV band, green from the 0.65-1.5 keV band, and blue from the 1.5-6 keV band. 



Fig. 5. — Color bars representing the possible colors obtainable from the combination of red, green and blue 
channels, referring to the mapped-color smoothed image in Fig. 4. The bars are labelled with corresponding 
values of HR1 and HR2. The tridimensional cube indicates the cuts in the values of blue counts we performed 
to produce the color bars (see text). 
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Fig. 6. — X-ray color-color diagram of HR1 vs. HR2. The blue grid indicates the regions occupied in the 
case of a power-law model, varying the intrinsic absorption Nh and the photon index T. The red grid 
represents the loci for a thermal model, varying the temperature kT and the metallicity Z. The filled square 
and triangle indicate the positions in the diagram of the regions described in §3. 
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Fig. 7. — Top: Mapped-color line-strength image of The Antennae. Red represents emission from the Fe-L 
line blend, green from the Mg line, and blue from the Si line (for details, see text). The 17 regions used for 
a preliminary subdivision of the hot ISM are marked in white. Bottom: 0.3-6 keV image of The Antennae, 
obtained removing the 3cr point sources. The 21 regions used for the spectral analysis of the hot ISM are 
marked in white. 



Fig. 8. — Chandra ACIS-S spectra for regions 1, 2, 3, 4a, 4b, and 5. Plotted for each region are the spectrum 
and the data-to-best-fit-modcl ratios. 



Fig. 9. — Chandra ACIS-S spectra for regions 6a, 6b, 7, 8a, 8b, and 9. Spectrum and data-to-best-fit-modcl 
ratios as in Figure 8. 



Fig. 10. — Chandra ACIS-S spectra for regions 10, 11, 12a, 12b, 13, and 14. Spectrum and data-to-best-fit- 
model ratios as in Figure 8. 



Fig. 11. — Chandra ACIS-S spectra for regions 15, 16, and 17. Spectrum and data-to-best-fit-model ratios 
as in Figure 8. 
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Fig. 12. — Distribution of the A% 2 values as a function of the temperature in four of our regions chosen for 
spectral analysis. The black dashed line in the upper-left panel represents the A% 2 distribution as a function 
of the lower temperature in a two-temperature best-fit model. The black dotted line, the red dotted line, 
and the green dotted line represent the 68%, 90% and 99% confidence levels, respectively. 
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Fig. 13. — Distribution of the A\ 2 values as a function of the Neon abundance in four of our regions chosen 
for spectral analysis. The black dotted line, the red dotted line, and the green dotted line represents the 
68%, 90% and 99% confidence levels, respectively. 
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Fig. 14. — Distribution of the A\ 2 values as a function of the Magnesium abundance in four of our regions 
chosen for spectral analysis. The black dotted line, the red dotted line, and the green dotted line represent 
the 68%, 90% and 99% confidence levels, respectively. 
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Fig. 15. — Distribution of the A\ 2 values as a function of the Silicon abundance in four of our regions chosen 
for spectral analysis. The black dotted line, the red dotted line, and the green dotted line represent the 68%, 
90% and 99% confidence levels, respectively. 
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Fig. 16. — Distribution of the A\ 2 values as a function of the Iron abundance in four of our regions chosen 
for spectral analysis. The black dotted line, the red dotted line, and the green dotted line represent the 68%, 
90% and 99% confidence levels, respectively. 
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Fig. 17. — Correlation between the measured EW and the abundance measured from the spectral fits for 
Mg (left) and Si (right). 
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Fig. 18. — Measured EW of the Mg line vs. measured EW of the Si line for the hot ISM regions of the 
spectral analysis. Right: the green lines represent the expected EWs as a function of the temperature 
(growing with growing values of EWsf, for details, see text) for a plasma in ionization equilibrium with 
Z = jZq, Z = Zq and Z = 5Zq. Right: the red lines represent the expected EWs as a function of the 
shock velocity (growing with growing values of EWsi] for details, see text) for a shock-heated plasma with 
Z = jZq, Z = Zq and Z = 5Zq 
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Tabic 1. Log of Chandra Observations 



Instrument 


Obs ID 


Date 


Grating 


Performed duration (ks) 


ACIS-S 


315 


December 1, 1999 


None 


75.53 


ACIS-S 


3040 


December 29, 2001 


None 


63.76 


ACIS-S 


3043 


April 18, 2002 


None 


60.81 


ACIS-S 


3042 


May 31, 2002 


None 


67.27 


ACIS-S 


3044 


July 10, 2002 


None 


36.49 


ACIS-S 


3718 


July 13, 2002 


None 


34.72 


ACIS-S 


3041 


November 22, 2002 


None 


72.91 
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Table 2. Net counts in the regions selected for the spectral analysis of the Diffuse Emission of The 

Antennae. 



R pen on 


CtSQ.3- 


0.65keV 


PTiSn ax. i t; iu ,-,!/" 

*-1"jU.dc>— l.ofce V 


.0 — D.[)KC V 


^Li^U.o — o-U/ce V 


1 


205. 


.3 


± 


16.4 


692.6 ±27.2 


71.4 ± 14.0 


969.1 ±34.6 


2 


162. 


.1 


± 


13.4 


565.8 ±23.9 


55.8 ±9.4 


783.9 ±29.2 


3 


218. 


.2 


± 


17.8 


754.3 ±28.8 


101.3 ± 17.0 


1074.1 ±37.9 


4a 


182. 


.9 


± 


14.5 


1114.8 ±33.7 


104.3 ± 12.5 


1402.0 ±38.7 


4b 


374. 


.0 


± 


20.4 


2556.0 ±51.0 


176.8 ± 16.0 


3107.0 ±57.2 


5 


176. 


.4 


± 


13.8 


876.5 ±29.6 


115.0 ± 11.9 


1167.9 ±35.0 


6a 


374. 


A 


± 


20.4 


1733.3 ±42.0 


231.4 ± 17.5 


2339.1 ±49.8 


6b 


520. 


.2 


± 


24.3 


3324.3 ±58.0 


286.7 ±20.7 


4131.7 ±65.8 


7 


115. 


.8 


± 


13.7 


1236.6 ±35.8 


379.1 ±22.7 


1731.3 ±44.4 


8a 


407. 


.9 


± 


20.8 


3772.4 ±61.7 


600.8 ±25.5 


4781.9 ±70.0 


8b 


313 


,6 


± 


18.3 


3128.4 ±56.0 


335.3 ± 19.4 


3777.3 ±62.1 


9 


714. 


.4 


± 


29.2 


5094.6 ±74.1 


703.3 ±30.9 


6510.2 ±82.3 


10 


583. 


.5 


± 


28.4 


2495.0 ±51.9 


271.1 ±26.7 


3349.8 ±65.0 


11 


216. 


.8 


± 


16.0 


1341.1 ±37.0 


174.2 ± 15.7 


1732.1 ±43.2 


12a 


338. 


,6 


± 


18.8 


1308.2 ±36.2 


81.7 ± 10.6 


1728.3 ±42.4 


12b 


153. 


.5 


± 


12.7 


662.1 ±25.9 


37.7 ±7.2 


853.5 ±29.6 


13 


119. 


.4 


± 


11.3 


561.3 ±23.9 


70.5 ±9.2 


751.0 ±28.0 


14 


53. 


.2 


± 


7.6 


633.3 ±25.1 


68.8 ±8.8 


755.1 ±28.0 


15 


131 


.6 


± 


11.8 


1031.3 ±32.1 


163.9 ± 13.3 


1326.7 ±36.6 


16 


482. 


3 


± 


24.3 


2530.8 ±51.0 


245.6 ±20.8 


3258.4 ±60.1 


17 


189. 


.7 


± 


18.6 


1296.3 ±37.4 


306.4 ±24.7 


1792.6 ±48.6 



Table 3. Best-Fit Parameters (lex errors for one interesting parameter) for the Diffuse Emission of The Antennae. 



rvcgion 


A/ 2 1 fin f 

X l ao J 


N H 

(X10 20 cm- 2 ) 


fcT 
(koV) 




Zm 9 

(xZ M g,Q) 


(xZsi.o) 


(X-ZFe,©) 


p 


IVIOQ 


1 


43.5/38 


< 2.68 


0.61+°°? 

— 0.03 


< 0.28 


0.23+" !J 

" — 0. IS 


< 0.19 


0.20+"™ 

-0.O4 




Al 


2 


34.3/25 


< 0.79 


0.62+"'™ 

— .04 


0.26 + „"o« 
— .26 


0.28+"^ 

" u — 0. 19 


0.08+" -f\t 

w ^ w — .08 


0.22 ± 0.04 




Al 


3 


57.8/48 


< 1.19 


0.54 ± 0.04 


0.39+1!" 

— 0.30 


0.47+^i 

— 0.23 


0.38+!!^ 

— 0.34 


0.23+ 1 ' 1 '! 

— 0.04 




Al 


4a 


51.1/44 


1.42 + ?J2 

— 1 .42 


0.63 ± 0.02 


0.29+"^ 

— 0.29 


0.80+''^ 

^" — 0.30 


0.27+!!' 22 

— 0.18 


0.38+''T 7 

— 0.07 




Al 


4b 


88.2/63 




0.20 ± 0.04 
0.60 ± 0.02 


- ^r»4-0 92 

1.631^.58 


„,-,4-0 97 


., „„4-0 75 

!-59-o.45 


°-96to 3 ii 




A2 


5 


31.4/35 




0.30 ± 0.02 


8.02+. 2b 8 * 


17.28ll!S 


-11.05 


2 76+ u - 7a 


1.88* 2 ' 


Bl* 2 ' 


6a 


74.6/61 


74+ i BU 

u - '^-0.74 


0.66 ± 0.02 


14+ U1U 

U- 1 -0.14 


0.28tS;iS 


24+" 14 

UlZ -0.12 


0.24 ± 0.03 


-0 61+"' 22 
u.o±_ 17 


Bl 


6b 


95.5/80 


1 47+116 


0.61 ± 0.01 


66 +U ' 21 

U.DO_ 19 


85+°- 1S 


70+' ) I7 


42+ IUM 




Al 


7 


86.9/65 


1J - "-4.05 


56+° 04 

U - OD -0.05 


2.68+^ 3 2 


2.93+-f 3 


o 97 + 3.31 
°- z '-1.68 


59+ iiu 


1.88 <2) 


B l(2) 


8a 


102.5/90 


8.6«+i:8il 


61+° 01 


±.OO_ 3 g 


0.87+.-- 


0.86+.-- 


n ro+U.OU 
0.t>8_ os 


1.88 (2 > 


Bl* 2 ' 


8b 


82.5/71 


24.61 + i^ 


0.33 ± 0.02 


±0 '^-6.39 


8-6itS:Si 


17.581-3°/ 


8.371";- 


3 R1+ 11 - 22 
°' D1 -0.09 


Bl 


9 


120.6/108 


n on + 1-18 
^■"^-1.00 


0.62 ± 0.01 


, on+0.33 
1 - za -0.29 


o.90i^; 22 


0.93± ; 22 


0.57 ± 0.07 


1.88 <2) 


B1 (2) 


10 


128.2/93 


< 0.31 


0.63 ± 0.02 


0.33 ± 0.16 


31+ uil 


„ Hi if; — 

0.57+^-;° 


0.27 ± 0.02 




Al 


11 


57.2/54 


1.88_ 1 88 


0.62 ± 0.02 


64+"- 4 ' 1 


1 43+ u;Ji 


1 02+"- 44 

1 - uz -0.35 


0.38_ 06 




Al 


12a 


55.0/46 


< 0.34 


0.60 ± 0.02 


15+ U1U 


0.24 ± 0.12 


°- 08 -0.08 


0.26 ± 0.03 




Al 


12b 


24.5/28 


< 0.1 . 


0.61 ± 0.02 


< 0.12 


20+°- iS 


< 0.14 


26+ " 5 




Al 


13 


32.2/22 


u - JZ -0.32 


58+" 04 


32 + u - 74 


°-°°-0.38 


62 + ' ) d7 


n 35+U .04) 
u.oo_ 08 


1.88 (2) 


Bl* 2 ' 


14 


28.8/21 


3 r>7+ 5 - 25 


0.37/pm0.03 


l-92-? 2 9 2 


- L -- L1 -0.80 


i.93± 2 ; 3 « 


0.65_ 18 


1.88 <2) 


B1 (2) 


15 


44.9/38 


o qc+4.08 


0.59 ± 0.03 


-1 ocj+1-24 


1 41+iM 


0.66tS;»» 


66+° 23 
u - uo -0 .18 


1.88 <2) 


B1 (2) 


16 


67.0/78 


0.76+i;^ 


0.62 ± 0.01 


31+ U ly 


54+°- ib ' 


59+" 17 


33+ lu)4 




Al 


17 


90.5/73 


3-08+ 2 ;^ 


62+' )u;i 


1 60+ iob 


2 QA+' 2 ' 2 " 


1 39+ 114 


°-8°-'o.2i 


1.88 (2) 


Bl* 2 ' 


(1) XSPEC 


Models: 


Al — wabs(wabs(vapec))\ A2 


— wabs(wabs(vapec+vapec))\ 


Bl = 


wabs(wabs(vapec+powerlaw)); 


B2 = 



wabs(wabs(vapec-l-vapec-f-powerlaw)). 
Photon index frozen. 



Table 4. Best-Fit Parameters (Z a and Zp e determined as described in the text and la errors computed for all interesting parameters) for 

the Diffuse Emission of The Antennae. 



Region # 



N H 

(XIO 20 cm-') 



kT 
(keV) 



(xZ Ne 



(xZm 9 , ( 



Z Si 
(xZ s . 



Z Fe 



1 


< 8.60 




< 0.75 




< 0.51 


o^otS™ 




2 


< 4.95 


„ -4-0 Ml 

o.62^:is 


„~4-ll N4 

o.2o+^;i* 


rrTTiTTn 

0-^035 


^"TrTTTTT*! 

O-19-o. 19 


„ 4-11 1 !i 

O-20+o.; 2 




3 


< 6.73 


-i-o rro 

u.o^_ 14 


; — -i- 1 oi 

0.48+ ; 48 


: -1-0 

U.OD_ 54 


: _|_| U| | 

o.5i+ =y 


m — rn 

U - ZD -0.11 




4a 


- — .~+8 io 

1A2 142 






„ „„4-() 76 




~ .,4-014 

44^ 
u -**-0.13 




4b 


97+ 4 ' 22 

u-y -0.97 


0. 20 (unconstr.) 
> 54 


1 or + 1.46 
1 " 5o -0.90 


1 44+3-93 


1 21 + 113 

1 -' i - L -0.80 


°- 8t> -0.34 




5 


n Q1+10.S9 
U - Mi -0.91 


o 30+ U07 


8 86+ 12 - 47 

O.OD_ 4 24 


19 r)7+ 2u 3-43 
LJA " -14.08 


iU9 -°' S -109.82 


2 81+ 8V - 4a 


1.88 (1) 


6a 


0.74+^ 


0.66I 07 


0-16+^16 


0.28+£« 


n nc + 0.47 

0.26 T 26 


0.27 ± 0.09 


-0.61 + ^ 


6b 


-1.47 


0.61 ± 0.04 


0.T2i-. 44 


94+"- 4u 

u ' a ^-0.37 


0.80_ 38 


43+"- 12 




7 


, „ --+11.69 
IO. i < -9.39 


56+ 42 


9 QK + 3-8& 
z - JO -1.83 


i - MO -1.19 


2 74 +3.u2 

z - '*-1.77 


46+"- 4;i 

U -^ D _0.18 


1.88* 4 ' 


8a 


8.66_ 4 81 


61+" 04 

u,u -0.05 


4 - &U -0.87 


U - OZ -0.53 


n R1+0.78 
°- 61 -0.44 


70+ u :il 

u - ' u _o.l9 


1.88( 4 ) 


8b 


24 61 + ll ti " 


ri r,o + U.2H 
u -°°-0.05 


1 5 57+°° 


5.95™ 


12-44^2.04 


tl - ii -4.35 


3 61+"- B4 

• 3 - Di -0.70 


9 




62+ " 4 


1 26+°' Jy 

i -^°-0.67 




0.84±«;S 


59+ 11 - 20 
u.oy_ 14 


1.88^) 


10 


< 2.04 


63+ " 4 


33+ 041 

U.OO_ Q 3 3 


41+°- 3i 


70+ u - 5(i 


25+ u u7 




11 




0.62 ± 0.06 


U - Oa -0.59 


, „-+0.94 
-0.75 


U - 8 '-0.63 


44+"- 17 




12a 


< 2.30 


60+ " 6 


11+ 048 


35+°- au 

u ' oo -0.32 


20 +u b0 


27+ u uu 

' -0.07 




12b 


< 5.29 


0.61 ± 0.06 


< 0.38 


14+ U(i5 
U.i"4_ 14 


< 0.82 


n oa +U.V2 

°- 26 -o.n 




13 


°-32 + 3 3 l 2 4 




2.57±»;# 


1 33+ 1Bls 


1 C4+13B.71 
- L - u ^-1.54 


47+^-88 

u ' -0.30 


l.SS* 1 ' 


14 


°' u ' -3.07 






unconstr. 


unconstr. 


unconstr. 


1.88 (1) 


15 


O qr + lU.92 


n 5Q+ U UU 


l-09i?; 9 




37+1-61 


72 + 1 2 3 
u - ' z -0.34 


l.SS* 1 ' 


16 


76 +4 " v 


62+' 1 " 4 


U -"-0.32 


0.52± U .f 2 


56+' 144 

U.OD_ 34 


34+' 14 ' 8 

u - -0.07 




17 


3 08+ 1026 


°- 62 -0.34 


96 +l4 V4 


2.91™ 


1 16+ 4 ' J4 

1 -- LU -1.16 


84 +3 - 25 

u '° -0.32 


l.SS* 1 ' 



to 



^ Photon index frozen. 
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Tabic 5. Best-Fit Parameters (68% errors) for the Integrated Emission of The Antennae (Comparison 

with ASCA). 



X 2 /dof 


(xlO 21 cm" 2 ) 


kT Z 
(keV) (xZ ) 


1043.5/371 




0.59 ± 0.03 0.18±8;S| 1.52±g;5| 
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